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Abstract 

In this paper, the iterative method developed by Daftardar-Gcjji and 
Jafari (DJ method) is employed for analytic treatment of Laplace equa- 
■ tion with Dirichlet and Neumann boundary conditions. The method is 

f^*) | demonstrated by several physical models of Laplace equation. The ob- 

lO ■ taincd results show that the present approach is highly accurate and 

| requires reduced amount of calculations compared with the existing iter- 

ative methods. 



Introduction 



- T I ■ 

^ | The Daftardar-Jafari method (DJM) developed in 2006 has been extensively 

■ used by many researchers for the treatment of linear and nonlinear ordinary 

and partial differential equations of integer and fractional order [IHZ]- An 
excellent modification of DJM can be found in [8] for solving different mod- 
els of linear and nonlinaer Klein-Gordon equations. The method converges 
to the exact solution if it exists through successive approximations. For 
concrete problems, a few number of approximations can be used for nu- 
merical purposes with high degree of accuracy. The DJM does not require 
any restrictive assumptions for nonlinear terms as required by some existing 
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techniques. The aim of this work is to effectively employ DJM to obtain ex- 
act solutions for different models of Laplace equation. While the variational 
iteration method [9HTT] requires the determination of Lagrange multiplier 
in its computational algorithm, DJM is independent of any such require- 
ments. Moreover, unlike the Adomian decomposition method [12H20] . where 
the calculation of the tedious Adomian polynomials is needed to deal with 
nonlienar terms, DJM handles linear and nonlinear terms in a simple and 
straightforward manner without any additional requirements. 

As mentioned before, we aim to obtain exact solutions of four models 
of Laplace equation, two with Dirichlet boundary conditions and two with 
Neumann boundary conditions. We also aim to establish that DJM is pow- 
erful, promising and efficient in handling nonlinear engineering problems. In 
what follows, we give a brief review of DJM. 



2 The DJ Method 



(2.2) u(x) = J2 Ul (H 



Consider the following general functional equation 
(2.1) u(x) = f(x) + N(u(x)) 

where N is a nonlinear operator from a Banach space B — > B and / is a 
known function, x = (x\,X2,--- ,x n ). We are looking for a solution u of 
(|2.ip having the series form 

oo 

_{x) 

n=0 

The nonlinear operator N can be decomposed as 

oo oo i i—1 

(2.3) N(J2 «i) = N(u ) + j> - N(£ Uj )}. 

i=0 i=l j=0 j=0 

From equations (|2.2p and (|2.3p . equation (|2.ip is equivalent to 

oo oo i i—1 

(2.4) i> = / + N ^ + E^Ek - N (E u ^ 

i=0 i=l j=0 j=0 

We define the recurrence relation 

{uo = f 
Ul = JV(uo) 
U m+ i = N(u H h Um) - N(u H h Um-x) 

Then 

(2.6) (ui H h u m ) = N(u H h Um) 



M. YASEEN ET.AL 



3 



and 

oo oo 

(2.7) Y,u l = f + N(J2^i) 

i=0 i=0 

The fc— term approximate solution of (|2.ip is given by u = uq+u\ + - ■ -+Uk~i 
In what follows, we apply DJM to four physical models to demonstrate the 
strength of the method and to establish exact solutions of these models. 

Example 1. Consider the second order laplace equation, 

(2.8) u xx -u yy = 0, < x, y < ir 
with boundary conditions 

, . u(x, 0) = sinh(x), u(x, n) = — sinh(x) 

u(0, y)=0, u(tt, y) = sinh(7r) cos(y) 

The equation ()2.8|) is equivalent to the following integral equation 

y y 

u = sinh(x) + yg(x) - [J (u xx (x, y))dydy 

o o 

where g(x) = u y (x,0). 

y y 

Set uq = sinh(x) + yg(x) and N(u) = - J f (u xx (x,y))dydy. 

o o 

Following the algorithm (|2.5h . the successive approximations are 



m = N{u Q ) = -^y 2 sinh(x) - ^y 3 g"(x) 

u 2 = N(u + «i) - N(u ) = ±y 4 smh(x) + ±y 5 g iv {x) 

u 3 = N(u + m + u 2 ) - N(u + ui) = -^y 6 smh(x) - ^y 7 g m {x) 



Combining the above results obtained for the components yields 

u(x, y) = smh(x)+yg(x)-^y 2 smh(x)-^y 3 g"(x) + ^y' 1 smh(x) + ^y 5 g m (x)+- ■ 

Using the inhomogeneous boundary condition, u(7r,y) = sinh(7r) cos(y) and 
the Taylor expansion of cos(y), we obtain 

111 y 2 y 4 

sinh(vr)+y5((7r)--?/ 2 sinh(7r)--yV'(7r) + — y 4 sinh(vr)+- • • = sinh(vr)(l- — + — 

Equating the coefficients of like terms on both sides gives 

g^) =g"(7r)=g™(n) = ■■■=() 

This means that g(x) = 0. Consequently, the solution in series form is given 
by 

u(x,y) = sinh(x)(l - — + — ) 
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and in closed form 

u(x, y) = sinh(x) cos(y) 
which is also the exact solution. 

Example 2. Consider the second order laplace equation, 

(2.10) u xx -u yy = 0, < x, y < 7r 

with boundary conditions 

/ 2 n -v u(x,0) = 0, u(a;,7r) = 

u(0,y) = sin(y), u(7r, y) = cosh(7r) sin(y) 

The equation (I2.10P is equivalent to the following integral equation 

y y 

u(x, y) = yg(x) - [J (u xx (x, y))dydy 
o o 

y y 

where g(x) = u y (x,0). Set u = yg{x) and N(u) = - f f (u xx (x,y))dydy. 

o o 

Following the algorithm (12. 5h . the successive approximations are 

ui = N(uq) = -±y 3 g"(x) 

u 2 = N(u + u 1 )-N(u ) = ±y 5 g iv (x) 

u 3 = N(u + m + u 2 ) - N(u + ui) = -jjy 7 g m (x) 



Combining the above results, we obtain 

u(x,y) = yg(x) - j/g"{x) + -^g iv {x) - ^y 7 g vi (x) ■ ■ ■ 

Using the inhomogeneous boundary condition, u(ir, y) = cosh(7r) sin(y) and 
the Taylor expansion of sin(y), we obtain 

y9(*)~yV(K)+^y 5 9™(K)~y 7 9 vi (K)+- ■ ■ = cosh(7r)(y-|^+|^- ■ ■ ). 
Equating the coefficients of like terms on both sides gives 
5 (tt) = /(Tr) = g iv (n) = ■ ■ ■ = coslnV) 

This means that 

g(x) = cosh(x) 
Consequently, the solution in series form is given by 

u(x, y) = cosh(x)(y - — + — ) 

and in closed form 

u(x, y) = cosh(x) sin(y) 
which is also the exact solution. 
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Example 3. Consider the second order laplace equation, 
(2.12) u xx - u yy = 0, < x, y < it 

with boundary conditions 

, , u y (x,0)=0, u y (x,n) = 2cos(2x) sinh(27r) 

1 ' u x (0,y) = 0, u x (7r,y) = 

The equation f)2. 12|) is equivalent to the following integral equation 

y y 

u = g{x) - (u xx (x,y))dydy 



o o 



y y 

where g(x) = u(x, 0). Set uq = g(x) and N(u) = — J J (u xx (x, y))dydy. 

o o 

Following the algorithm (|2.5p . the successive approximations are 



Ul = N{u ) = -^y 2 g"(x) 

u 2 = N(u + Ul ) - N(uo) = iy 4 g iv (x) 

u 3 = N{u + Ul + u 2 ) - N(u + ni) = -±y 6 g m (x) 



Combining the above results, we obtain 

(2.14) u(x,y) = g(x) - ^y 2 g"(x) + ^y 4 g m (x) - ^y 6 g m (x) + ■■■ 

Using the inhomogeneous boundary condition, u y (x, ir) = 2 cos(2x) sinh(2 
and the Taylor expansion of sinh(2-7r), we obtain 

1 1 87T 3 S2-7T 5 

-Trg"(x) + -^g n '(x) - -n 5 g vi (x) + • • • = 2 cos(2x)(27r + — + — + • • 

Equating the coefficients of like terms on both sides gives 

g"(x) = -4cos(2x) 

which gives 

g{x) = cos(2x) + Co 
Consequently, the solution in series form is given by 

4y 2 16m 4 

u(x,y) = cos(2x)(l + + + • • • ) + C 
and in closed form 

u(x, y) = cos(2x) cosh(2y) + Co 
which is also the exact solution. 
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Example 4. Consider the second order laplace equation, 
(2.15) u xx -u yy = 0, < x, y < it 

with boundary conditions 



(2.16) 



u y (x,0) = cos(x), u y (x,ir) = cosh(7r) cos(x) 



u x (0,y) = 0, u x (n,y) = 

The equation f|2. 15|) is equivalent to the following integral equation 

y y 

u = g{x) + y cos(x) - J j (u xx (x, y))dydy 
o o 

y y 

where g{x) =u(x,0). Set uq = g{x)+y cos(x) and N(u) = — j j \u xx (x,y))dydy. 

o o 

Following the algorithm (|2.5p . the successive approximations are 

ui = N(uq) = -^y 2 g"(x) + ^y 3 cos(x) 

u 2 = N(u + m) - N(uq) = \y A g iv {x) + ±y 5 cos(x) 

u 3 = N(u + ui+ u 2 ) - N(u + ui) = -hy 6 g m (x) + ly 7 cos(x) 



Combining the above results, we obtain 

(2.17) u(x,y) = cos(x) S mh(y) + g(x)-^y 2 g"(x) + ^y A g iv (x)-^y (i g vi (x) 

Using the inhomogeneous boundary condition, u y (x,ir) = cosh(-7r) cos(x) , 
we obtain 

cos(x) cosh(-7r) — irg"(x) H — -Tv 3 g iv (x) 77T 5 g m (x) + • • • = cosh(-7r) cos(x) 

3!! 6! 

Equating the coefficients of like terms on both sides gives 
g"(x)=g lv (x)=g m (x) = --- = 

which gives 

g(x) = C 

Consequently, the solution in series form is given by 

y3 y5 

u(x, y) = cos(x)(y + — + — H ) + C 

and in closed form 

u(x, y) = cos(x) sinh(y) + Co 
which is also the exact solution. 
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